QCD corrections to T production via color-octet states at the Tevatron and LHC 



O 

(N 

Or 
<D 

o 



43 ■ 

Oh 1 
i ' 

43 . 



> ■ 
os : 

m . 
oo ■ 
cn ; 

Os : 

o : 
o : 



X 



Bin Gong 1,2,3 , Jian-Xiong Wang 1,3 and Hong-Fei Zhang 1 ' 3 
Institute of High Energy Physics, CAS, P.O. Box 918(4), Beijing, 100049, China. 
Institute of Theoretical Physics, CAS, P.O. Box 2735, Beijing, 100190, China. 
Theoretical Physics Center for Science Facilities, CAS, Beijing, 100049, China. 

(Dated: September 21, 2010) 

The NLO QCD corrections to T production via S-wave color-octet states Y ^Sq 8 ', ^J 8 '] at the 
Tevatron and LHC is calculated. The K factors of total cross section (ratio of NLO to LO) are 1.313 
and 1.379 for T[^ 8) ] and T[^ 8) ] at the Tevatron, while at the LHC they are 1.044 and 1.182, 
respectively. By fitting the experimental data from the DO, the matrix elements for S-wave color- 
octet states are obtained. And new predictions for T production are presented. The prediction for 
the polarization of inclusive T contains large uncertainty rising from the polarization of T from feed- 
down of \b- To further clarify the situation, new measurements on the production and polarization 
for direct T are expected. 
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I. INTRODUCTION 

For heavy quarkonium production and decay, a naive 
perturbative QCD and nonrclativistic factorization treat- 
ment is applied straightforwardly. It is called color- 
singlct mechanism (CSM). To describe the huge dis- 
crepancy of the high-p t J/ip production between the 
theoretical prediction based on CSM and the experi- 
mental measurement at Tevatron, a color-octet mecha- 
nism [1] was proposed based on the non-relativistic QCD 
(NRQCD) [2]. In the application, J /if) or T related pro- 
ductions or decays are very good places for two reasons, 
theoretically charm and bottom quarks are thought to 
be heavy enough, so that charmonium and bottomonium 
can be treated within the NRQCD framework, experi- 
mentally there is a very clear signal to detect J/ip and T. 
The key point is that the color-octet mechanism depends 
on nonperturbative universal NRQCD matrix elements, 
which is obtained by fitting the data. Therefore various 
efforts have been made to confirm this mechanism, or to 
fix the magnitudes of the universal NRQCD matrix ele- 
ments. Although it seems to show qualitative agreements 
with experimental data, there are certain difficulties. A 
review of the situation could be found in Refs. [3, 4]. 

To explain the experimental measurements [5, 6] of 
J/ip production at the B factories, a series of calcula- 
tions [7, 8] in the CSM reveal that the next-to-leading 
order (NLO) QCD corrections give the main contribution 
to the related processes. Together with the relativistic 
correction [9], it seems that most experimental data for 
J/ip production at the B factories could be understood. 
Recent studies show that the NLO QCD correction also 
plays an important role in J/ip production at RHIC [10] 
and the hadroproduction of \c [H]- For the J/ip photo- 
production, the pt and z distributions can be described 
by the NLO calculations in CSM [12] by choosing a small 
renormalization scale, but recent NLO calculations in 
CSM [13] show that the pt distributions of the production 
and polarization for J/ip can not be well described when 
choosing a proper renormalization scale. Although the 



complete calculation at NLO in COM [14] can account for 
the experimental measurements on the pt distribution, it 
cannot extend to J/ip polarization case. To further study 
the heavy quarkonium production mechanism, there are 
many other efforts performed, such as NLO QCD correc- 
tion to J/tp production associated with photon [15], QED 
contributions in J/tp hadroproduction [16], inclusive J/ip 
production from T decay [17], double heavy quarkonium 
hadronproduction [18], and NLO QCD correction to J/ip 
production from Z decay [19]. 

For the polarized heavy quarkonium hadroproduction, 
the leading order (LO) NRQCD prediction gives a sizable 
transverse polarization for J / ip production at high p t at 
Tevatron while the experimental measurement [20] gives 
slight longitudinal polarized result. The discrepancy was 
also found in T production. In a recent paper [21], the 
measurement on polarization of T production at Teva- 
tron is presented and the NRQCD prediction [22] is not 
coincide with it. Within the NRQCD framework, higher 
order correction is thought to be an important way to- 
wards the solution of such puzzles. Recently, NLO QCD 
corrections to J / ip and T hadroproduction have been cal- 
culated [23-27], and the results show that the NLO QCD 
corrections give significant enhancement to both total 
cross section and momentum distribution for the color- 
singlet channel. This would reduce the contribution of 
color-octet channel in the production. Also, it is found in 
Ref. [24] that the polarizations for J/ip and T hadropro- 
duction via color-singlet channel would change drasti- 
cally from transverse polarization dominant at LO into 
longitudinal polarization dominant in the whole range of 
the transverse momentum p t at NLO. It seems that these 
results open a door to the solution of the problem. But 
things are not always going as expected. The NLO QCD 
corrections to the J/ip production via S-wave color-octet 
states were studied in our previous work [28]. It was 
found that the effect of NLO QCD correction is small 
and the discrepancy holds on. For the color-singlet part, 
the partial next-to-next-to- leading order (NNLO) calcu- 
lations for T and J / ip hadroproduction show that the un- 



2 



certainty from higher order QCD correction [38] is quite 
large, therefore no definite conclusion can be made. As 
we know, the contribution from the color-octet states is 
smaller in the T production than that in J/ip production, 
thus things may be different. In this paper, we present 
our calculation on NLO QCD corrections to T hadropro- 
duction via S-wave color-octet states. New matrix ele- 
ments are fitted and new prediction for the polarization 
status is presented. 

This paper is organized as follows. In Sec. II, we give 
the LO cross section for the process. The calculation of 
NLO QCD corrections are described in Sec. III. In Sec. 
IV, we present the formula in final integration to obtain 
the transverse momentum distribution of T production. 
Sec. V. is devoted to the description about the calcula- 
tion of T polarization. The treatment of J/ip is discussed 
in Sec. VI. The numerical results are presented in Sec. 

VII, while the summary and discussion are given in Sec. 

VIII. In the Appendix, several details of the calculation 
are presented. 

II. THE LO CROSS SECTION 

According to the NRQCD factorization formalism, the 
inclusive cross section for direct T production in hadron- 
hadron collision is expressed as 

a[pp^T + X] = J dx 1 dx 2 G i/p G j/p 

i,j,k,n 

X a[i + j^(bb) n + k](0%), (1) 
I 

di 

da B (gq^r[ 3 s[ S) ]q) 

di 

d<j B (gg^r[ 3 s[ 8) ]g) 
di 

da B (qq^T[V 8) ]g) 

di 

dfr B ( gg ->TpS< 8) ]g) 
di 

d& B (gg^?[ty 8) ]g) 
di 



by introducing three dimensionless kinematic variables: 

. _ (pi +P2) 2 f _ (pi - Pz) 2 „ _ (pi - Pi) 2 
8 Am 2 4 ™6 ' U ~ 4m b ' 

(3) 



where p is either a proton or an anti-proton, the indices 
k run over all the partonic species and n denotes 
the color, spin and angular momentum states of the in- 
termediate bb pair. The short-distance contribution a 
can be perturbatively calculated order by order in a s . 
The hadronic matrix elements (O^) are related to the 
hadronization from the state (bb) n into T which are fully 
governed by the non-perturbative QCD effects. In the 
following, a represents the corresponding partonic cross 
section. 

At LO, there are three partonic processes: 



9(pi) +9(P2) -+ T[^ 8) ,^ 8) ](p 3 ) +g(p A ), (Ll) 
g(pi) + q(P2) -+ r[ 1 si 8 \ 3 S[ 8) ](p 3 ) + q(p 4 ), (L2) 
q(Pi) + q(P2) -+ T[^ 8 \%p>](p 3 ) +g(p 4 ). (L3) 



where q represents a sum over all possible light quarks or 
anti-quarks: u, d, s, c, u, d, s and c. In our calculation 
of T production, we take charm quark as light quark as 
an approximation. Typical Feynman diagrams for these 
three processes are shown in Fig. 1. And the partonic 
differential cross sections in n = 4 — 2e dimension for LO 
processes can be obtained as 



(2) 

I 

and the reasonable approximation My = 2m& is taken. 
Our LO results are consistent with those in Ref. [29]. 
The LO total cross section is obtained by convoluting 



n'aUOZCSiMi-l) 2 + (g - l) 2 ][4t 2 - tu + 4u 2 ] 
324m5s 2 (s- l)Hu 

-K 2 al(Oj(%))[(s-l) 2 + (u-l) 2 ][As 2 -su + Au 2 ] 



864mgs 3 (i- l) 2 u 

Tr 2 a 3 s (Oj( 3 Si))[(s 2 - l) 2 + (i 2 - l) 2 + (w 2 - l) 2 - 6stu - 2] [19 - 27(si + tit + us 



+ 0(e), 

5 - 9,7(.$i -A- in -4- ii.*}] 

+ 0(e), 



1152m 5 b s 2 {t- l) 2 (u - 1) 2 (J - l) 2 

57r 2 a3(Oj(^ ))[i 2 + « 



+ 0(e), 
+ 0(e), 

4- Sf.ii] \s 4 4- f 4 4- u 4 4- 1 1 

0(e), 



216mp 3 (s- l) 2 
-^ 2 o? s {Oj{\))[s 2 +u 

576mgs 2 i(i- l) 2 
^ 2 al{Oj{\))[s 2 t 2 + s 2 u 2 + i 2 u 2 + siu][s 4 + i 4 + u 4 + l] 



256m 5 b sHu(t - l) 2 (u - l) 2 (s - l) 2 




FIG. 1: Typical Feynman diagrams for LO processes, a) Feynman diagrams for process (LI); b) Feynman diagrams for processes 
(L2) and (L3). Diagrams in groups (ai), (02), (61) and (62) are absent for the Sq 8 ' state. 
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the partonic cross section with the parton distribution 
function (PDF) in the proton: 

t< 8 > + x] = J2 J + 3 -> T(8) + fc ] ( 4 ) 

x Gi/ P (xi, fj,f)Gj/ p (x2,fJ>f)dxidx2, 

where denotes certain color-octet T[ 1 S'q 8 ' ) ] or 

, Hf is the factorization scale. 

III. THE NLO CROSS SECTION 

The NLO contributions can be written as a sum of 
two parts: first is the virtual corrections which arise from 
loop diagrams, the other is the real corrections caused by 
radiation of a real gluon, or a gluon splitting into a light 
quark-antiquark pair, or a light (anti)quark splitting into 
a light (anti) quark and a gluon. 

A. Virtual corrections 



where 7^ is the Euler's constant, /3q = ^-Ca — ^TpUf is 
the one- loop coefficient of the QCD beta function and n/ 
is the number of active quark flavors. We have four light 
quarks u, d, s and c in our calculation, so nf—A. The 
color factors are given by 7> = 1/2, Cf — A/3,Ca = 3 
and fj, r is the renormalization scale. 



There are 267 (for the state) and 413 (for the 3 sf ) 
state) NLO diagrams for process (LI), including counter- 
term diagrams, while for both processes (L2) and (L3), 
there are 49 (for the !?^ 8) state) and 111 (for the 3 s{ 8) 
state) NLO diagrams. Part of the Feynman diagrams for 
these processes are shown in Fig. 2. The diagrams in 
which a virtual gluon line connects the quark pair pos- 
sess Coulomb singularities, which can be isolated and 
attributed into renormalization of the bb wave function. 



For each process, by summing over contributions from 
all diagrams, the virtual correction to the differential 
cross section can be expressed as 



There exist ultraviolet (UV), infrared (IR) and 
Coulomb singularities in the calculation of the virtual 
corrections. UV divergences from self-energy and trian- 
gle diagrams are canceled by introducing renormaliza- 
tion. Here we adopt the renormalization scheme used in 
Ref. [30]. The renormalization constants Z m , Z2, Z21 and 
Z3 which correspond to bottom quark mass m&, bottom- 
field tpb, light quark field tp q and gluon field A a ^ are defined 
in the on-mass-shell (OS) scheme while Z g for the QCD 
gauge coupling constant a s is defined in the modified- 
minimal-subtraction (MS) scheme: 
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where MS -, is the amplitude of process (Lj) at LO, and 
Mj^ .j is the renormalized amplitude of corresponding pro- 
cess at NLO. MX.-i is UV and Coulomb finite, but it still 
contains IR divergences. And the total cross section of 
virtual contribution could be written as 



[pp -> T< 8 » + X] = J2 Jv V [i + 3^ T< 8 ' + k] (7) 
x Gi/ p (xi,(j,f)Gj/ p (x2,(Xf)dxidx2, 




FIG. 2: Typical one-loop diagrams, a) Feynman diagrams for gq — > T[jSq ]g and qq —¥ T^Sq ]g; a + b) Feynman diagrams 



for gq 



T[ 3 s[ S) ]q and qq 



]<?; c) Feynman diagrams for — > Y ['ISq 8 '] g; c + d) Feynman diagrams for gg 



Counter-term diagrams, together with corresponding loop diagrams, are not shown here. 



B. Real corrections 

There are eight processes involved in the real correc- 
tions: 



99 -> 


T[^ 8) 


3 c(8 
' a l 


]gg, 


(Rl) 


9Q -> 


T[^ 8) 


3 c(8 


]gq, 


(R2) 


qq^ 


T[^ 8) 


3 c(8 
) D l 


]gg, 


(R3) 


99 -> 
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3„(8 
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]qq, 


(R4) 


qq^ 


T[^ 8) 


3„(8 
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]qq, 


(R5) 


qq^ 


T[^ 8) 


3„(8 
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(R6) 


qq -> 


T[^ 8) 


3„(8 
) D l 


]qq, 


(R7) 


qq 1 -> 


T[^ 8) 


3„(8 
' °1 


}qq', 


(R8) 



where q, q' denote light quarks (anti-quarks) with dif- 
ferent flavors. Feynman diagrams for these processes 
are shown in Fig. 3. We have neglected the contribu- 
tions from the another two processes, gg —> T^bb and 
qq —> T^bb, which are IR finite and small. Phase space 
integrations of above eight processes generate IR singu- 
larities, which are either soft or collinear and can be con- 
veniently isolated by slicing the phase space into differ- 
ent regions. We use the two-cutoff phase space slicing 
method [31], which introduces two small cutoffs to de- 
compose the phase space into three parts. Then the real 
cross section can be written as 



It is easy to observe that different parts of IR singularities 
from one real process may be factorized and each part 
should be added into the cross sections of different LO 
processes. This is the reason why we have to calculate 
the NLO corrections to the three LO processes together. 

1. soft 

Soft singularities arise from real gluon emission. Thus 
only real processes (Rl), (R2) and (R3) contain soft sin- 
gularities, corresponding to the three LO processes. One 
should notice that, unlike color-singlet case, the soft sin- 
gularities caused by emitting a soft gluon from the quark 
pair in the S-wave color octet exists and we find that 
the factorized matrix element is the same as the case of 
emitting a soft gluon from a gluon. 

Suppose ps is the momentum of the emitted gluon. If 
we define the Mandelstam invariants as = {pi + Pj) 2 
and tij — (j>i — pj) 2 , the soft region is denned in term of 
the energy of p§ in the p\ + P2 rest frame by < E5 < 
8 S \Jsi2/2. For each of the three real processes, a s from 
the soft regions is calculated analytically under the soft 
approximation. 

Following the similar factorization procedure as ap- 
plied in the calculation of color-singlet case [25] , the ma- 
trix elements for a certain real process (Ri) in the soft 
region can be written as 

|M [Ri] | 2 | soft c-4™^ J2 ~ P \ Pk Mg p 

fki\ ■ P5)(Pk ■ P5) lM 

(9) 



Oi) (02) (a 3 ) (04) (a 5 ) (61) (6 2 ) 




(<fe) (ds) ('/.') (</:,) (A) (e) (7,i (/ 2 ) 

FIG. 3: Feynman diagrams for real correction processes, a) Rl (T [^ 8) ] ); a + b) Rl (T[^ 8) ]); c) R2~R4 (T[^ 8) ]); c + d) 
R2~R4 (T [%f >] )■ e) R5~R8 (T [^ 8) ] ); e+f) R5~R8 (T [%f ] ] ). Rl (T [&< 8) ] ) denotes process gg -> T [&< 8) ] 33, Rl (T [^5[ 8) ] ) 
denotes process gg —> T ["^Sj 8 '] gg, and so on. 



with 



2. hard collinear 



<] 



T°(j)M 



b 1 -b 1 



t r 



T a (fc)M^; 



(10) 



where Mr£T is the color connected Born matrix cl- 
ement for LO processes (Li). If the emitting parton 
j is an initial state quark or a final state antiquark, 
T a (j) = T b a ih .. For an initial state antiquark or a final 

state quark T°(j) = —T^. b i . If the emitting parton j is 

a gluon or the color-octet state, T a (j) = if a bjb;- And 
the corresponding parton level differential cross section 
can be expressed as 



da [Ri] = 



with 



r(l-e) fiirtf 



2ttT(1 - 2e) V S12 



j.k=i 



[Li] M 



(11) 



(12) 



The factor P k is universal for all three real processes, 
and is given in Appendix. (A). Sometimes doh^ may be 
written in a more compact form as 



da- 



[Li] [L« 



(13) 



where Ct -, is a constant. This is always true if the LO 
process (Lj) contain only one independent color factor 
in the matrix element. But for processes with two or 
more than two independent color factors, there seems no 
sure reason for it to be or not to be true. Of course, 
no matter Eq. (13) is true or not, we can always obtain 
dot , through Eq. (12). Most processes involved in this 
calculation have more than one independent color factors, 
and they are listed in Appendix. (B 1). 



The hard collinear regions of the phase space are those 
where any invariant (s^ or iy) becomes smaller in mag- 
nitude than S c si2- It is treated according to whether 
the singularities are from initial or final state emitting or 
splitting in the origin. 

a. final state collinear For real processes (Rl) ~ 
(R6), which contain final state collinear singularities, the 
final state collinear region is defined by < S45 < 5 c Si2- 
Again following the similar factorization procedure de- 
scribed in Rcf [31], the parton level cross section in the 
hard final state collinear region can be expressed as 



~HC\n 1 -Br T n 



r(l-e) flirtf 



2itT(1 - 2e) V S12 



.4 



HC 



(14) 



For a certain real process (Rj), (L-) is the corresponding 
LO process it factorizes into. And the coefficient A^ c 
are listed in Table. I, with 



•99 _ 



j\9~^99 



Aq~^ Q9 

j\9-+qq 
j[9-^qq 



and 



•99 _ 



N (11/6 + 
N [67/18 
C F (3/2 ■+ 
C F [7/2 - 
-n//3 
n//3 (In S, 



6' 



2lnS' s ) 

-7r 2 /3-ln 2 S' 3 -lnS c (ll/6 + 2lnS' 8 )] 
2\nS' s ) 

ir 2 /3-\n 2 5' s -ln5 c (3/2 + 21n^)] 



5/3) , 



(15) 



A' 12 



Sl2 + S45 



Mi 



(16) 
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Thus the total cross section for real correction processes 
in hard final state collincar region can be written as: 



with 



HC 



where 
B HC (g) 



E [af c [i + j^rw + k 1 + k 2 } 

i,j,ki,k 2 

xG i/p (x 1 , Hf)G j/p (x 2 , fif)dx 1 dx 2 
E J a B [i + j^T^+k]B HC (k) 

xGj/p^i, Hf)G j/p (x2, Hf)dx 1 dx 2 , (17) 
a s r(l-e) fAntf 



Af c = 2A sc (g -> gg) 

A S 2 C = A sc (q^qg) + A sc (g^gg) 

Al c = 2A sc (q^qg), (20) 



and 



B HC (q) 



27rr(l - 2e) V s 12 

^ + A 9 ^ 99 + A 3 -* qq 

27rr(l-2e) ^ s 12 J \ \ e 



(18) 



+ 4™ 









gg -»■ T< 8 ) 55 




_j_^9^99 i ^9~>99 


gg T^pg 


gq -+ T< 8 >g 


—A 9 + A 9 


55 -> T< 8 >gg 


5<? -»• TWj 


f ^ - ^ - 

_j^9^11 _|_ ^9^<3<3 


gg ->■ T (8) ff3 


gg^T< 8 'g 


_ ^9-^99 _|_ ^9->99 


gg — > T^'gg 


gg^T< 8 > 5 




gg^T< 8 'g'g' 


gg-^T< 8 ' 5 





TABLE I: The hard final state collinear factors for real cor- 
rection processes and the corresponding LO processes. 



b. initial state collinear Almost all real processes, 
except process (R6), contain hard initial state collinear 
singularities. These singularities are partly absorbed into 
the redefinition of the parton distribution function (PDF) 
of the concerned hadrons (usually it is called as the mass 
factorization [32]). Here we adopt the scale dependent 
PDF using the MS convention given in Ref [31]. 



G b/p (x,fi f ) = G b/p (x) 



a s r(l - e) / 4irfj% 



27rr(l-2e) V fij 



x / — P w {z)G v/p {x/z). (19) 

J x z 

The second term is sometimes referred as the mass fac- 
torization counter-term. There is still something remain- 
ing after the cancellation, which can be expressed in two 
terms. The first one, which only exists in the real pro- 
cesses with final state gluon, can be expressed as 



af c [R,] = a»\Li] 



a s T(l - e) f Anul 



27rr(l - 2e) V M/ 



A 



sc 



A sc (q^qg) = - e C F [3/2 + 2ln(5 s )] 
A sc (g^gg) = J [2N\a8 s + {UN - 2n f )/6] .(21) 

The corresponding hadronic total cross section is 

E Jtrii+j^r^+k+g] 

xG i / p (xi,fj,f)G j / p (x2,fif)dx 1 dx2 

E )^ +3 -+ T(8) + fc HS 5C (^) + B sc (j)} 



HC 



x G i/p (x i , (if ) Gj/p (x 2 , (i f ) dxi dx 2 , 



(22) 



with 

B sc (g) = 
B sc (q) = 



a s r(f — e) / 47T/X 



.27rr(l - 2e) V M/ 
a s r(l-e) 
2ttT(1 -2e) \ (if 



A sc {g^gg) 

A SC {q^qg). 



The other term is obtained by summing up the remaining 
contributions from all the real correction processes. It 
can be written as 



T (8) + k 



(23) 

a s r(l-e) /4^\ e 



27rr(l -2e) \ s 



Gi/ P (xi, (if)Gj/p(x2, (if) + {xi «-» x 2 ) 



dxidx 2 , 



with 



1-6,5,. 



■ dy, 



G c/p (x,nf) = E / —G C '/p(x/y,(if)Pcc'(y) , 

(24) 

and 



P ij (y) = P ij (y)ln(5 c ^-y S -^)-P( j (y). (25) 

The n-dimensional unregulated (y < 1) splitting func- 
tions Pij(y,e) has been written as Pij(y,e) = Pij(y) + 
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ePr(y) with 

P qq (y) 
P' qq {y) 
P 9q (y) 
KM 
P 99 (y) 
KM 
P q9 (y) 
p q9 (y) 



c F 



i + y 2 



1-2/ 
Cf(i-»). 

1 + (1 - y) 2 



2iV 



2/ 



+ +2/(1-2/) 



1-2/ 

[2/ 2 + (1 - 2/) 
-2/(1-2/)- 



(26) 



C. Cross section of all NLO contributions 

The hard noncollinear part a HC is IR finite and can be 
numerically computed using the standard Monte-Carlo 
integration techniques. Now the real cross section can be 
expressed as 



a =a +a f + a i + cr add + a 



And we have 



N LO B , V _,_ R 



(27) 
(28) 



IV. TRANSVERSE MOMENTUM 
DISTRIBUTION 

To obtain the transverse momentum p t distribution 
of T, a similar transformation for integration variables 
(da^di — > Jdp t dy) which we introduced in our previous 
work [25] is applied. Therefore we have 



with 



J Jdx^yd/pixx, Hf)G j/p {x 2 , fJ. f )^-, (29) 



Pi =zi — (1,0,0,1), p 2 =x 2 — (1,0,0,-1), 



m t = ^M T +Pt, P3 = (m t cosh2/,p t ,0,™ t sinhy) 
2m t 



Xi\ 



y/S' 

4xix 2 pt 
2x\ — x t e y ' 

_ 2r + x t e y 
~ 2-xte-y' 



X 2 



2r + x\ x t e~ v 
2x\ — x t e y 



(30) 



where y/S is the center-of-mass energy of pp(p) at Teva- 
tron or LHC, 771,4 is the invariant mass of all the final 
state particles except T, and y and p t are the rapidity 
and transverse momentum of T in the laboratory frame 
respectively. 



POLARIZATION 



The polarization parameter a is defined as: 



da T /dp t - 2da L /dp t 
d<r T /dp t + 2da L /dp t ' 



(31) 



It represents the measurement of T polarization as func- 
tion of T transverse momentum p t when calculated at 
each point in p t distribution. To evaluate a(p t ), the po- 
larization of T must be explicitly retained in the calcula- 
tion. The partonic differential cross section for a polar- 
ized T is expressed as: 



p i -e(X)p j -e*(X), (32) 



where A = T U T 2 ,L. e(Ti), e(T 2 ), e(L) are the two trans- 
verse and longitudinal polarization vectors of T respec- 
tively, and the polarizations of all the other particles are 
summed over in n-dimension. One can find that a and 
&2 j tire finite when the virtual corrections and real cor- 
rections are properly handled as aforementioned. There- 
fore there is no difference in the differential cross section 
da\/dt whether the polarization of T is summed over in 4 
or 77 dimensions. Thus we can just treat the polarization 
vectors of T in 4-dimension, and also the spin average 
factor goes back to 4-dimcnsion. The gauge invariance 
is explicitly checked by replacing the gluon polarization 
vector into its 4-momentum in the final numerical calcu- 
lation. 



VI. TREATMENT OF J/tp 

The production mechanism of J/ip at Tevatron and 
LHC is much similar to that of T except that, color-octet 
states contribute much more in J / ip production according 
to the experimental data and LO theoretical predictions. 
The results of above calculation can also be applied to 
the case of J/ip by doing the substitutions: 

777b ^ fn c 
R S (0) T i? s (0) J ^ (33) 

77/ = 4 «-> 77/ = 3 

Note that in J/ip production, charm quark is no longer 
treated as light quark. 
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VII. NUMERICAL RESULT 

In our numerical computations, the CTEQ6L1 and 
CTEQ6M PDFs [33], and the corresponding fitted value 
a s (M z ) = 0.130 and a s (M z ) = 0.118 are used for LO 
and NLO calculations respectively. The bottom quark 
mass is set as 4.75 GcV. 

The choice of the renormalization scale p r and factor- 
ization scale /j,f is an important issue in the calculations, 
and it causes uncertainties. We choose p = p r = pf — 
\J (2m&) 2 + p\ as our default choice. And the center-of- 
mass energies are chosen as 1.96 TeV at Tevatron and 14 
TeV at LHC. 

At First, different values of the two cutoffs, S s and S c , 
are used to to check the independence of the final re- 
sults on the cutoffs and the invariance is observed within 
the error tolerance. Then the two phase space cutoffs 
are fixed as S s = 10~ 3 and S c = S s /50 in the following 
calculations. 

It is known that the QCD perturbative expansion is 
not good in the regions of small transverse momentum or 
large rapidity of T. Therefore, the results are restricted 
in the region p t > 3. For the rapidity cut, \yy\ < 1-8 
is chosen at the Tevatron, the same cut condition as the 
experiments [34], and at the LHC, it is chosen to be \y\ < 
3. 

To fix the NRQCD matrix elements for color-octet 
states of T(IS'), the DO data [34] is used, and the fitting 
starts from Eq.(4) of Ref. [35] where the contributions 
from spin-singlet states r]i,(nS) and hb(nS) are not in- 
cluded. And we have to take a few approximations in 
our fitting procedure: 

• For the S-wave color-singlet part, only the direct 
color-singlet T(IS') and feed-down from T(2S) are 
considered, while other contributions have been ne- 
glected. The contribution from the feed-down of 
T(2S) can be included to the direct T(IS') pro- 
duction by multiplying a factor of Br[T(2S) — > 
T(IS') + X] x (of(2S))/(Oj (15)), which results 
in a factor of 1.127 after a short calculation with 
PDG data [36]. And the results for direct T(1S) 
of color-singlet contribution are extracted from our 
previous work [25]. 

• The contributions from P-wave color-singlet states 
Xb.j(nP) are estimated by multiplying a decay frac- 

tion F ^J } » + F Xb \J y where F^ 

and ^j(2P) can De obtained from an older sam- 
ple with the cuts p t > 8 and \yy\ < 0.4 [37]. As 
pointed out in Ref. [38], the fraction should not 
depend very strongly on p t according to Fig. 2 of 
Ref. [39]. Also, from Fig. 4 of Ref. [34] we can see 
it should not depend very strongly on the rapid- 
ity cut either. Thus F r Q S l = 27.1 ± 6.9 ± 4.4% 

and F T( * S „\ = 10.5 ± 4.4 ± 1.4% are taken in our 

calculation, which result in F \ L ss 37.6 ± 9.4%. 



• The contribution from P-wave color-octet states 
T[Py'} at NLO are still not available. As shown 

below, the NLO QCD corrections to T[!s^ 8) ] don't 
change the cross section very much. If we assume 
that the NLO QCD corrections to T^j 8 -*] are also 

small, we can mix it with T[ 1 5q 8 ' ) ] again, like what 
we have done at LO. Thus, the value of our fit- 
ted (Oj(So)) inc includes the contributions from 
T[ 3 pj 8) ] as well. 

With these approximations, the formula we used for the 
fitting of inclusive color matrix elements becomes 

da[T] inc = 1.127 x&p)!^)]^^)) 

+ Fli%d<r[T] inc + d ( T[(^) 8 ( 1 5 )](0 8 T (!So)) inc 
+ rfa[(66) 8 ( 3 5 1 )](0 8 T ( 3 ,9 1 )) inc , (34) 

and the NRQCD matrix elements for color-octet states 
(Oj) inc are determined as 

{OjCSo)).^ = (0.948 ± 0.444) x 10~ 2 GcV 3 
(OjCSi)).^ = (4.834 ±0.719) x 10~ 2 GcV 3 , (35) 

where only the uncertainty in F^^fl, has been consid- 

Xb ) 

ered. The fitting is shown in Fig. 4, together with our 
prediction for inclusive T production at the LHC. 

The direct fraction of direct T production can also 
obtained from Ref. [37] as Fj^ = 50.9 ± 12.2%. Thus 
we can use the formula 

F d T r (1S W[T] inc = daKbbU's^OjCSr)) 
+ da[(bb)sCSo)](Oj(%)) 
+ da[(66) 8 ( 3 ,9 1 )](0 8 T ( 3 5 1 )), (36) 

to fit the direct color-octet matrix elements. The matrix 
elements are obtained as 

(OjCSo)) = (0.630 ± 0.576) x 10~ 2 GeV 3 

(Oj Cfci)) = (3.900 ± 1.063) x 10~ 2 GeV 3 , (37) 

where the uncertainty comes only from Fj^ 1 ^. Again 
the value of our fitted (OJ^Sq)) includes the contribution 

from T^j 8 -*]. This fitting is shown in Fig. 5, together 
with our prediction for direct T production at the LHC. 
The band in the figure is obtained from the uncertainty 
of F T(1S) 

The dependence of the total cross section on the renor- 
malization scale n r and factorization scale p, / are shown 
in Fig. 6. It is obvious that the NLO QCD corrections 
make such dependence milder. We can also see that the 
NLO QCD corrections effect the cross section lesser at 
the LHC than at the Tevatron. The pt distributions of 
T production via S-wave color-octet states are presented 
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FIG. 6: Total cross section of T hadroproduction at LHC 
FIG. 4: Transverse momentum distribution of inclusive T pro- ( upper curves ) and Tevatron (lower curves) , as function of /x 
duction at Tevatron and LHC. The DO data is from Ref. [34] . ^ ^ = ^ = ^ and ^ = ^ {2 m h y+pl 









1 DO data 


1 






total LHC 

tota 

W8) 

Jo 

3 S (B) 

3 s< 


10 








_Q 








c 








a_ 2 

^ 10 

t- 

m 

X 




— + — 




CL -3 
T3 10 








b 








-4 

10 









6 8 10 12 14 16 IS 20 



P, (GeV) 

FIG. 5: Transverse momentum distribution of direct T pro- 
duction at Tevatron and LHC. The DO data is from Ref. [34]. 



in Figs. 7 and 8, where only slight changes appear when 
the NLO QCD corrections are included. 

T[ 1 5q 8 ^] produces unpolarized T, so it contributes to 
a = for both LO and NLO. The pt distributions of 
T polarization parameter a from Y[S' n p' ) ] are shown in 
Fig. 9 and there is slight change when the NLO correc- 
tions are taken into account. Our predictions for the 
polarization of direct T production are also presented in 
the figure as a "total" result. In Fig. 10, the polarization 
of inclusive T production at the Tevatron is shown. As 








5 10 15 20 25 50 55 40 45 50 

Pt(GeV) 

FIG. 7: Transverse momentum distribution of T production 
with fi r — [if = /in at the Tevatron. 



the polarization of T from the feed-down of Xb(nP) is not 
available yet, a huge band is obtained by verifying the po- 
larization of this part between -1 to 1. The experimental 
data from the DO is also shown in the same figure. We 
can see that, there is still some distance between the the- 
oretical prediction and experimental measurement, even 
with such a large band. 
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FIG. 8: Transverse momentum distribution of T production 
with fi r = fj,f ■ — fio at the LHC. 
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FIG. 9: Transverse momentum distribution of polarization 
parameter a for direct T production at the Tevatron (upper) 
and LHC (lower). 



VIII. SUMMARY AND DISCUSSION 

As a summary, in this work, we have calculated the 
NLO QCD corrections to T production via S-wave color- 
octet states Tf^Q 8 ^, 3 S[ 8S> ] at the Tevatron and LHC. 
With fi r = fif = /i0) the K factors of total cross section 

(ratio of NLO to LO) are 1.313 and 1.379 for T[^ 8) ] 

and T[% 1 (8) ] at Tevatron, while at LHC they are 1.044 
and 1.182 respectively. Unlike for the color-singlet case, 
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FIG. 10: Transverse momentum distribution of polarization 
parameter a for inclusive T production at the Tevatron. The 
DO data is from ref [21]. 



there are only slight changes to the transverse momen- 
tum distributions of T production and the T polariza- 
tion when the NLO QCD corrections are taken into ac- 
count. All the results imply that the perturbative QCD 
expansion quickly converges for T production via the S- 
wave color-octet states, in contrast with that via color- 
singlet, where the NLO contributions are too large to 
hint a good convergence at the NNLO. By fitting the 
experimental data from the DO at the Tevatron, the ma- 
trix elements for S-wave color-octet states are obtained. 
And new predictions for the pt distributions of the T 
production and polarization at the Tevatron and LHC 
are presented. The prediction for the polarization of in- 
clusive T contains large uncertainty rising from the po- 
larization of T from feed-down of Xb- Even with such a 
large uncertainty, there are still some distance between 
the prediction and experiment data. Also, the errors of 
fractions used in the fitting, F^ip)' ^Xbi^P) anc ^ ^cUr > 
are quite large and result large uncertainty in the matrix 
elements. New measurements on the production and also 
polarization for direct T are expected. This would make 
the matrix elements more precise, and get rid of the large 
uncertainty from \b- 

This work is supported by the National Natural Sci- 
ence Foundation of China (No. 10475083, 10979056 and 
10935012), by the Chinese Academy of Science under 
Project No. INFO-115-B01, and by the China Postdoc- 
toral Science foundation No. 20090460535. 



Appendix A: Calculation of the factor I-* 

If we write the n- momentum of soft gluon in the p\ +P2 
rest frame as 

p 5 = E 5 (l, . . . , sin6>i cos 9 2 , cos Q x ) , (Al) 

then P k is defined as 

J [Pj ■ P5)(Pk ■ P5) 
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with 
dS 



7T \ S 12 

x sin 



2e e 2 de 2 . 



dE 5 E\- 2t sin 1 - 2 ^! d6i 
(A3) 



Before the calculation of P k , define 0j as 0j — \pj \/Ej, 
which is the ratio of momentum to energy of particle i in 
the pi + pi rest frame, where 

Pi = fa= fa = 1, 3 = ^ = 0. (A4) 

s + 1 

Then we can write pj and pk as 
Pj = E^l,--- ,0j) 

p k = E k (l,- ■ ■ ,0 k sm6 jk ,0 k cos6j k ), (A5) 
where 6j k is the angel between j and k. Now we have 
T jk _ 1 - 3 0k cos0 jk 

1 — 1 E 1 A 

where 

.S\2j Jo 



If= — 



dE 5 Ei 



(A6) 



-l-2£ 



5.) W 



(A7) 



and 

P A k = fwi 1 - 2e 1 dB 1 r sm- 2f e 2 de 2 - — ^ — - 

Jo J a 1 - ft cos6»i 

x \ (A8) 

1 — 0u cos9 cos^i — 0k sinO sin#i cos# 2 

The way to calculate the integrals P A can be found in 
the appendix of Rcf. [31]. Now we come to the results. 
It's easily to obtain 

jii = /22 = /44 = Qj (AQ) 

and the others are listed below. 

1. I u with i = 2,3,4. Write the momenta of the 
particles as 

pi = ,1), 

P2 = S 2 (l,--- ,-1), 

P3 = S 3 (l, • • • ,/3sinfli3,/3cosfli 3 ), 
p 4 = S 4 (l,--- ,-sin0i 3 ,-cos0i 3 ), (A10) 
then we have 



rl2 
i A 



rl3 



rl4 



7T 

! 

e 



, -+ln 

1 — cos 6*13 ^ e 

-iln 2 J + 2Li 2 (£)-2Li 2 ( 

^ \ 1 — t 



ln 2 (l - t) 



(All) 



2tt 



(1 + cos 6*13)6 \1 — s 



1 + e 2 Li 2 



1-s 



which lead to 

rl2 ^ ,12 1 r-2e 

7 = --lEl A = ~-^ d s 



rl3 



7r e 
1 — /3cos0 13 



l.'i 



1 - eln 



2e 2 s 1 

1 9 

-ln 2 J + 2Li 2 (t)-2Li 2 
^ VI — t 



ln 2 (l-£) 



14 1+CQS6>13 r .14 
i = 1 E^A 



-^5 



(A12) 



2 "s 



1 - S 



1 + e 2 Li 2 



1 - s 



2. I 2t with z = 3,4. These two can be directly ob- 
tained from In with the substitution t o u. 



ln 2 (l-u) 



1 



--In s + 2Li 2 (u)-2Li 2 . 
2 '" \ 1 — u 



i 24 = -i<5r 2e/ ; 



1 - s 



1 + e 2 Li 2 



1 - s 



3. I 33 and I 34 . Write the momenta of the final state 
particles as 



Ps = E 3 (l,--- ,-0), 
p A = Ei(l,--- ,1), 



(A13) 



then 



r 33 



r34 
1 A 



I 33 -- 



2tt 



l-^ 2 

7T 



1 + e^lns 



l + 0\-e +lnS - e 



ln 2 s + 2Li 2 (l-s) 



I-/? 2 



/3 4 = _ 



7T 
1+0 



T r 3 3 _ A~2e 



1 + e-lns 



t r 34 



(A14) 



1 

2?' 



1 - e In s + ^ 



ln 2 s + 2Li 2 (l-s) 



Appendix B: Color Factors 

Here we present color factors for all the processes in- 
volved. Color indices for particle n are labeled as j n . 



1. LO processes 

The color factors listed here for LO processes have been 
orthogonalizcd and normalized. 
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gg — > T^S^jg, three color factors in total: 
1 

^ - Tr l~T-' 4 T-' 3 T-' 2 T- 71 — T-' 4 T- 71 T-' 2 T- 73 



-Tr[r :,4 r : ' 2 r : ' i r : ' 3 — t 1i t 33 t 3i t 32 



75 

1 

7! 



gg T[S[ ]g, three color factors also: 
1 



3V2 



3^ 



Tr [(T 3i T 32 T 33 T 31 + T 3i T 3l T 33 T 32 

_ (rpil r pil r p02 r pj3 _|_ r p3A'J<33'J<32rpjl 
^Uj'3i r ph r pj3 r pjl _|_ J<3iJ<3l r p33 r p32 
_|_ (T<3i r p3l r [<32 r [<h _|_ J 1 ^ JT?3 yj2 

— 2(T J4 T J2 T J1 T 33 + T^T^T^T^ 

^\^p34 r p32 r pj:-i r pji _|_ r pji r pji r pj3 r p]2 

_|_ ^ r pji r ph r [ ) h r ph _|_ J>J4 J>J3 J>J2 J>jl 
_|_ ^['34.'J-'32'J'3irjlj 3 _|_ rpj 4 rpj 3 rpj 1 rpj 2 

• gq — > TT^q 8 ^, only one color factor: 
' [3(T^T^+T^T^) j4j2 -^ 2 ^ 3 ]. 



(Bl) 



2V15 1 

• .gg — > Y[ 3 S^]g, two color factors: 



(B2) 
(B3) 



^(T*T^« a ,-^(8T*T* +T»T*) jija . (B4) 
• gg — > T^g, almost same with gg — > T^q. 



2. Virtual correction processes 

In the amplitude of virtual correction processes, be- 
sides the same color factors as in the corresponding LO 
process, there are extra ones. As we have mentioned be- 
fore, virtual correction to the cross section is related to 
virtual amplitude as Eq.(6). Then the terms in propor- 
tion to these extra color factors will vanish and do not 
contribute to the final result as we have orthogonalized 
the color factors of LO processes. Thus no new color fac- 
tors in virtual correction processes need to be presented 
here. 



3. Real correction processes 

In order to present the color factors of real correction 
processes in a simplified form, we list here all indepen- 
dent color factors in a certain process. Actually in our 
calculation, they are orthogonalized and normalized too, 
which are too complicated to be listed here. 



• gg — > T\S^]gg, twelve color factors. The permu- 
tations of ji, j2, j3 and j4 contain 24 terms. Divide 
them into twelve groups with two terms in each 
group, and the twelve color factors can be expressed 
as 

Tr [T-* 5 (T a T 6 T c T d + T d r c r & T a )] (B5) 
where a,b,c,d are permutations of j\ , ji , j'3 and j\ . 

• gg — > Y^S^] gg, also twelve color factors. They 
can be expressed as 

rjlj. ^rparpbrpCrpd _ rpdrpCrpbrpCl -^ (B6) 

• gg —> T [^S'q 8 ''] gg, five independent color factors: 

d hj3k (T jl T k ) j4h , d 3233k (T k T 31 ) jih , 
d jlj3k (T j2 T k ) j4j5 , d jlj3k (T k T j2 ) jij5 , (B7) 
Q(T j3 T j2 T jl -T jl T j2 T j3 ) j4j5 +if jlj2j3 5 j4j5 . 

• gg — > Y^S^jgg, seven independent color factors. 
One is d Jlj2j3 Sj 4 j 5 while the others can be expressed 
as 

(T a T b T c ) Mj57 (B8) 
where a, b, c are permutations of ji, j'2 and j'3. 

• gq — > T^gq and gg — > T^qg, similar to — > 
T^gg. 

• gg — > T[S , Q 8 ^]gg, two color factors: 
Tll h 5nn+Tlt h 5 hh , Tj 3 ^ +T^8 hh . (B9) 

• gg — > Y^j^jgg, four color factors: 

T njJoiji^ T j4j 1 <>32hi Tjlji&iijsi T hjjnh- ( B l°) 

• gg — > Y^gg, similar to gg — > T^'qq. 

• gg — » TfS'Q 8 ^] g'g', only one color factor: 

^ T njJoiOi + T j4j 1 $32h) - 2 ( T j2h S 3d 5 + T ujJ]23i)- 

(Bll) 



gg -4- T[^ 8 ^] g'g', three color factor: 
T 33 d T 33 - ■ 

°- L jdl U 3235 1 jih°323l- 



(B12) 



gg' 



' -> T( 8 )gg', similar to gg -> T^q'q'. 
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